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The importance of accounting for an ambipolar electrostatic potential or a non-uniform density 
profile for modelling of inductively coupled discharges is demonstrated. A drastic enhancement is 
observed of the power transfer into plasma for low-collisional, low-pressure, non-local discharges 
with non-uniform electron density profiles under the condition of bounce resonance. This enhanced 
plasma heating is attributed to the increase of the number of resonant electrons, for which the 
bounce frequency inside the potential well is equal to the rf field frequency. 
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INTRODUCTION 

Low pressure radio-frequency (rf) inductive discharges 
have been extensively used over the past decade as 
sources of inductively coupled plasmas (ICP) in the 
plasma aided material processing industry, semiconduc- 
tor manufacturing, and lighting [jj, |2| ■ For very low pres- 
sures, i.e. in the milliTorr region the ICP discharges ex- 
hibit a strong non-local behavior and a number of pe- 
culiar physical effects typical for warm plasmas, like an 
anomalous skinpenetration and a resonant wave-particle 
interaction [3, |j] . The study of these effects leads to fur- 
ther optimization of the ICP sources and can result in 
improvement of the characteristics of plasma-based de- 



An interesting effect that can lead to enhanced heat- 
ing for bounded, low-pressure plasmas is the possible 
bounce resonance between the frequency u> of the driv- 
ing rf field and the frequency of the bounce motion of 
the plasma electrons confined in the potential well by an 
ambipolar potential 4>{x) and the sheath electric fields 
near the discharge walls 0, 0, 0, H H E3, El- Most 
earlier theoretical and numerical studies on this subject 
assumed for simplicity a uniform plasma density over the 
discharge length, and the absence of an ambipolar poten- 
tial, i.e. allowed the electrons to bounce inside a poten- 
tial that is flat inside the plasma and infinite at the walls 
H El El E, El Although these suppositions can result 
in a fairly good description of the plasma behavior under 
non-resonant conditions, the discharge parameters under 
resonant conditions can be greatly altered by account- 
ing for the presence of the ambipolar potential, which - 
it should be stressed - always exists in real discharges. 
It is a very well known result of the quasilinear theory, 
that for low-collisional discharges the plasma heating es- 
sentially depends on the so-called "resonant electrons," 
or electrons with velocities equal to the phase velocities 
v ~ u/k of the plane waves constituting the rf field (Lan- 
dau damping) j^, 0]. For bounded plasmas and for an 



electron to be resonant, the above condition transforms 
into the requirement that the rf field frequency must be 
equal to, or be an integer multiple of the bounce electron 
frequency ui — nfi&. But the electron bounce frequency is 
very sensitive to the actual shape of the ambipolar elec- 
trostatic potential </>(x), especially for low-energy elec- 
trons. Accounting for the electrostatic potential can lead 
the plasma electrons into the resonant region even if they 
were not there in the absence of the potential. This can 
result into a drastic enhancement of the plasma heating 
and other related phenomena 

In this article we present the results of a full, self- 
consistent numerical modelling of low pressure ICP dis- 
charges under the bounce resonance condition and show 
the pronounced influence of the presence of the electro- 
static ambipolar potential on the plasma parameters un- 
der resonant conditions. 



BASIC EQUATIONS 

The model assumes a one-dimensional, slab geome- 
try, inductively coupled discharge of a plasma bounded 
on both sides by parallel walls with a gap length L. 
The walls carry fixed currents flowing in opposite direc- 
tions, produced by an external radio frequency source. 
The induced solenoidal rf electric field E y is directed 
along the walls and the static ambipolar electric field 
E x = —dtp/dx of the ambipolar potential 4>{x) is directed 
towards the discharge walls, keeping electrons confined 
and the plasma quasineutral, i.e. n e (x) = rii(x). In 
the present treatment of high density discharge plasmas 
(n e ~ I0 8 — 10 12 cm -3 ) the sheath width is neglected, be- 
cause it is of the order of a few hundreds of microns, much 
less than the discharge dimension L. Furthermore, it is 
assumed that the plasma electrons experience specular 
reflection: a) from the discharge walls when they have 



total energy e — mv /2 



(x) larger than the elec- 



tron potential energy —e<f>(x w ) at the walls, ±x w , and 
b) from the geometrical location of the turning points 
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x ±{ £ ), where — e<j>(x±) = e. The above 1-D scheme can 
also be a good approximation for a cylindrical ICP dis- 
charge, if the rf field penetration depth S into the plasma 
is less than the plasma cylinder radius R . 

In order to describe the discharge self-consistently, one 
needs to determine the rf electric field profile E y (x), the 
electron energy distribution function (EEDF) fo(z), and 
the ambipolar potential 4>(x). The detailed description of 
all the needed formalism is given in 0] . A short account 
of the formalism is given below. 



Calculation of the EEDF 

For low-pressure discharges, where the energy relax- 
ation length is large compared with the plasma width 
and the energy relaxation time is large compared with 
the rf period, the electron velocity distribution func- 
tion (EVDF) can be represented as a sum of the main 
isotropic part / = fo(e) (EEDF) that is a function of only 
the total energy e and of a small alternating anisotropic 
part/ifov,*), f = fo{e) + h{x,v,t) The 
Boltzmann equation for the electron velocity distribution 
function reads 



dfx 
dt 



dh 
dx 



e d(p dfx 
m dx dv x 



eE y (x,t) djfp + h) 
m dv y 
= St(h + f ), 



(1) 



where E y (x,t) is the nonstationary rf electric field, and 
St(f) is the collision integral. After applying the stan- 
dard quasilinear theory, Eq.QJ splits into two equations 
0, a linear one for fx 



dfx 



dfx ed4>dfi eE y (x,t) df 



dt dx m dx dv x 
and a quasilinear one for /q 



in 



OVy 



eE y (x,t) dfx 



dv„ 



St(fo). 



St(fx), (2) 



(3) 



Here, the bar denotes space-time averaging over the 
phase space available to electrons with total energy e 
|l7l flil Il9j . We can represent as harmonic functions 
the rf electric field E y (x,t) = E y0 (x) exp(— iwt) and the 
anisotropic part of the EVDF f\ — fiQ exp(— iu>t), where 
lj is the discharge frequency. Using the Bhatnagar-Gross- 
Krook (BGK) approximation yj|, St(fx) = —vfx> and 
omitting the subscript in the amplitudes, Eq. can 
be rewritten as 



IWfx + v x-q^\ 



ev y E y (x) — 



-vfx 



(4) 



where v is the transport collision frequency, e x — 
mv^,/2 + (p(x) is the total energy along the a;-axis, and 
(p(x) — —e(j)(x) is the electron potential energy. 



Eq. Q can be effectively solved using a Fourier series 
expansion. Introducing the variable angle of the bounce 
motion 0] which is proportional to the time of flight 
of an electron from the left turning point to the current 
point 



I X, S X 



Trsgxx(v x ) f x dx 



T{e x ) J x _ \v x (e x )\' 



(5) 



where T is half of the bounce period of the electron mo- 
tion in the potential well ip(x) 



T b (e x ) 



x + dx 



x_ \vx{e x )Y 



(6) 



Eq. 10} simplifies to 



jfi + n b ^\ ea -v y eE y (6)^ = -vfi. (7) 

where tlbi^x) = ir/T(e x ) is the bounce frequency for the 
electron in the potential well. Making use of the Fourier 
series 



l(x,e x )= ^ g n exp(inO), 



9n 



2tt 



g(0, e x ) exp (—inn 9) d6 



Eq. J7J gives 

E y n (E-x) 

and 



E y (0) cos (n8)d0 



(8) 



(9) 



(10) 



f ls (x,e x ) = l/2(fx(v x >0) + fx(v x <0)) (11) 

dfo 
de 



-mvyVy (x,s x ) — , 



where 



e 
rn 



E - 



Eyn cos[n0(x) 



in^lt, — itu + v 



(12) 



Knowing the symmetrical part f\ s of the anisotropic 
contribution to the EVDF, one can average Eq. J2J ac- 
cording to 



Term(a;, v)(e) = / dxv(x,e)Term[x,v(x,e)],(13) 

J X- 

v(x,e) = y/2[e - ip(x)]/m. (14) 
and obtain the final equation for /o 

-±(D e + D-)^--^\vr e + V^fo= (15) 



E 



y*(w + e%) V^JM/ ( £ + £ *) - K f Q 
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Here, the bar denotes averaging according to Eq. (|13fl . 
and is the inelastic collision frequency. The coefficients 
V e i , D ee , V ee stem from the elastic and electron-electron 
collision integrals, respectively, and are given by flill2lj| 



2m 

Vei = ^rrwv 



M 



2wv„ 



dw^/wf , 



(16) 



(17) 



Note that the normalization factor in Ea. (|23|l appears 
due to the normalization of /o as 



n e (x) = / f (e)^E- tp(x)de. 

Jcp(x) 

We now use the Fourier series 

oo 

E y( x ) — E Ss cos(k s x), 



(24) 



(25) 



4 wv,, 



3 n 



(J dww 3/2 f + w 3 ' 2 J dwij, (18) 



47rA e 



(19) 



where s is an integer, k s — (2s + 1)tt/L. Substituting 
Eq. J211) into Eq. yields 



-kj + —\ H„ = — 



21 ' 



(26) 



where w = mv 2 /2 is the electron kinetic energy, v ee is 
the Coulomb collision frequency, and A ee is the Coulomb 
logarithm. 

The energy diffusion coefficient responsible for the elec- 
tron heating is given by 



2 oo ~ e 



(20) 



£ - e, 



n— — oo 



^(e*) [n 6 (e x )n - + z/ 2 



Note that this expression for D £ (e) accounts for the 
bounce resonance Qb(sx)n — The dependance of elec- 
tron plasma heating on resonant electrons especially pro- 
nounced for the v << u, as in this case 



[Clb(e x )n — v] + 1/2 



ir5(Q b n — u>) 



(21) 



where S() is a Dirac delta function. It is worth to note 
that if L — > oo, the summation in l|21|l goes into integra- 
tion over corresponding wave vectors k n , and the bounce 
resonance condition Q b (e x )n — to transforms into the 
well-known wave-particle resonance condition for contin- 
uous wave spectrum kv = lj. 



Calculation of the rf electric field 



The transverse rf electric field E y is obtained from a 



single scalar equation 



d 2 E y | u 2 E 
dx 2 c 2 v 



\j(x)+I5(x)-IS(x-L)}, 



(22) 

where / is the wall current and j{x) is the induced elec- 
tron plasma current density that can be calculated know- 
ing the anisotropic part f\ s of the EVDF 



an 3 / 2 

3 = ~- 7^ I flsVyd^V. 



Att^/2 



(23) 



e n e \p „ „gen 

mi{2s + l)fi 6T s - 1 \{2s + l)n bT 



, (27) 



where n e is the plasma density at the discharge center, 
n bT = V T n/L, V T = y/2T/m, and we introduced the 
generalized plasma dielectric function 



/2(2s + l)7rf7 



bT 



E 



nVL b {e) - (2s + l)Q, bT £, ft 6 (e) 



where 



T(e) 



fo(e)de. 



(28) 



de, 



(29) 



In the limit of a uniform plasma, the generalized dielec- 
tric function coincides with the standard plasma dielec- 
tric function [15j. The coefficients Gi !n (s) are the tempo- 
ral Fourier transforms of cos(fc/:r) in the bounce motion 
of an electron in the potential well (md 2 x/dt 2 = edep/dx) 



Gi, n {e) = i 



cos[kix(r)] cos (J-^-^j dr 



(30) 



The Maxwell equation l|26l) together with the equation for 
the electron current lj2"?| and ptljl comprise the complete 
system for determining the profiles of the rf electric field. 



Calculation of the electrostatic potential 

The electrostatic potential is obtained using the 
quasineutrality condition 



n e (x) = rii(x) = / f (e)yje- tp(x)de, (31) 

Jip{x) 
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where n e (x) is the electron density profile and rii(x) is 
the ion density profile given by a set of fluid conservation 
equations for ion density and ion momentum 



drij ^ dinjUj) 
~dt 



dx 



and 



d(niUi) d(niUiUi 



at 



0:t 



n, d(j){x) 
"Mi dx 



(32) 



(33) 



where R is the ionization rate, Vi is the ion-neutral colli- 
sion frequency and rii,Ui, and Mi arc ion density, velocity, 
and mass, respectively. 

Eq. I|3f I) is solved in the form of a differential equation 
IH 



dt P T scr/ ^ d\n[m(x)} 

— = —l e (X) ; . 

dx e w dx 



where T| cr (a;) is the electron screening temperature 

-l 



TT\x) 



1 



Me)- 



de 



2n{x) J v{x) y/ £ - ip(x) 



(35) 



and the electrostatic ambipolar potential can be obtained 
by integration of Eq. 1)34( 1 . 

The above described self-consistent system of equa- 
tions was formulated in [15j, and implemented and com- 
pared with the experimental data in [22I l23j | . Although 
the simulation results of the latter articles were proven 
to be adequate, the method of the direct computation of 
Green functions, used there, seems to be impractical, be- 
cause of the excessively long computational time (about 
a day). To speed up the calculations (to about an hour), 
the Fast Fourier Transform [Eas. 1)26)1 -1)30 )1 ] was used in 
the present simulations. 



RESULTS AND DISCUSSION 

Collisionless heating is a very important channel of 
power transfer for bounded, warm, low-collisional plas- 
mas . For a semi-infinite plasma, the collisionless heat- 
ing essentially depends on the resonant electrons, i.e. 
electrons moving with the velocities equal to phase ve- 
locities of the components of the spectrum of the driv- 
ing rf field. Such electrons can effectively gain energy 
from a wave and, henceforth, the plasma can be effi- 
ciently heated under resonant conditions. For the case 
of bounded plasma, the condition for resonance heating 
transforms into the bounce resonance condition, i.e., the 
frequency w of the electromagnetic wave must coincide, 
or to be several times larger of the bounce frequency f2b 
of the electron bounce motion in the potential well. If the 
electron mean free path is larger than the discharge gap 
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FIG. 1: The electron bounce frequency Qt(£x) = n/Tb(e x ), 
for discharge of length L — 5 cm, as a function of the 
(34) electron energy [s x = mv x /2 — e<f>(x)] for different poten- 
tial wells, consisting of the reflecting walls and different am- 
bipolar potentials 4>{x). Solid line corresponds to a uniform 
plasma with <f>(x) = 0; dashed line - quadratic potential 



eV; and dotted line - quartic potential 

+ 



<p(x) = 5 x (2x/L 

(t>(x) = 5x (2x/L-l) 4 eV. Here, T b (e x ) = dx/v x ,e x (e x ) 
The cross-hatched box shows the resonance region. Arrows 
show electron energies in the resonance region. 



L, the resonant electrons with = lo accumulate veloc- 
ity changes in successive interactions with the rf electric 
field, which lead to very effective electron heating [jj. 

The importance of the resonant electrons for plasma 
heating under the bounce resonance condition, can be 
readily examined from expression 1)21)1 for the electron 
energy diffusion coefficient, where the term Qb(£x)n = w 
in the denominator clearly shows the dominant role of the 
resonant electrons for collisionless electron heating. The 
total power P deposited into plasma, per unit square of 
a side surface, is related to the electron energy diffusion 
coefficient D £ (e) [l5| . 



P = 



De[£) d JM de . (36) 
10 de 



The presence of an ambipolar electrostatic potential 
can greatly affect the electron heating due to two reasons: 

1) the ambipolar potential confines low energy elec- 
trons to the center of the discharge plasma and these 
electrons cannot reach the region of the strong field near 
the walls, and 

2) the number of resonant electrons is generally larger 
for a nonuniform plasma than for a uniform plasma due 
to influence of the electrostatic ambipolar potential on 
the bounce frequency. 

If the ambipolar potential is accounted for, then for low 
energy electrons, the distance between turning points is 
smaller than the distance between walls L. This results 
to an increase of their bounce frequencies compared to 
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FIG. 2: The plasma resistivity as a function of the distance 
between the walls for a uniform plasma (without any ambipo- 
lar potential) and a nonuniform plasma with the quadratic 
and quartic potentials and a given Maxwellian EEDF. Dis- 
charge parameters are: electron temperature T e = 5 eV, 
peak electron density at the center of the discharge n e — 
5 x 10 11 cm" 3 , rf field frequency cj/2tt = 13.56 MHz, and 
electron transport frequency v = 10 7 s _1 . The lines corre- 
spond to the same cases as in Fig. 1. 
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the uniform plasma case. Therefore, if electrons with 
low energy were far from the bounce resonance in a uni- 
form plasma, they can approach the resonance region in 
nonuniform plasma. Fig. ^ shows the dependance of the 
electron bounce frequency f2j, (e x ) on the electron energy 
e for different potential wells, consisting of the reflect- 
ing walls and different ambipolar potentials <j>(x). Here, 
Qb(£x) = 2n/Tb(e x ), where Tt is half of the bounce period 
of the electron motion in the potential well ip(x) given by 
Eq. ®. The width of the resonance is given by Ea. pijl 
and is proportional to v. The population of resonant elec- 
trons consists of all electrons corresponding to the inter- 
val of bounce frequencies Q,b(e x )n £ [ui — v, ui + v] where 
for most practical cases n = 1 - only the first resonance is 
important. From Fig.^it is apparent that the number of 
resonance electrons increases if the ambipolar potential 
is accounted for. For example, all electrons confined in 
the quadratic potential have the same bounce frequency 
and are all resonant, if the bounce frequency is equal to 
the discharge frequency. The results of the calculations 
described below show that increase of the number of the 
resonant electrons has a much more profound effect on 
the discharge plasma heating than the mere confinement 
of them in the region of low rf field. 



FIG. 3: The plasma resistivity as a function of the driving 
frequency normalized by (a) the inverse of the transit time 
through the skin layer Vr/5max and (b) by the bounce fre- 
quency Vrn/L, for different electron densities and for the 
fixed discharge parameters: L = 50 cm, T e — 5 eV and 
v = 10 6 s -1 . Here, Vt is the thermal velocity and <5 max is 
the rf field plasma penetration depth for the maximal value 
of the resistivity at a given density. 



Non-self-consistent simulations with a given 
Maxwellian EEDF and given ambipolar potentials 

Most of earlier works on low-pressure ICP discharges 
where reported assuming a Maxwellian EEDF, uniform 
plasma and accounted for the sheath electric field by 
bouncing electrons off the discharge walls |8|, llfll llll 1 1 2| . 

In order to explicitly show the importance of account- 
ing for ambipolar potential on collisionless electron heat- 
ing, we performed numerical simulations using a given 
Maxwellian EEDF for uniform and nonuniform plas- 
mas (with and without an ambipolar potential). Specif- 
ically, we obtained results for the dependence on the 
plasma length of the plasma resistivity or the real part 
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of impedance Z = Att/c x E(0)/B(Q), where E(0) and 
-B(O) are the electric and magnetic field at the wall, re- 
spectively |3- The surface impedance is related to the 
power deposition according to: 

P = -I 2 ReZ, (37) 

where / is the amplitude of the current. 

The results are presented in Fig. |2 It is clearly seen 
that the presence of the ambipolar potential enhances sig- 
nificantly the resistivity of the plasma under the bounce 
resonance condition, compared to the case of the ab- 
sence of a potential. The most profound change in re- 
sistivity is observed for the quadratic potential. In this 
latter case all trapped electrons have the same bounce 
frequency, and thus all of them can be resonant. The 
obtained results explicitly show that neglecting the am- 
bipolar potential, as is often done for simplicity, can lead 
to large discrepancies (more than 100 percent), espe- 
cially for conditions close to the bounce resonance. For 
large gaps L ^S> 2tt6 = 5 cm (where the skin depth is 
S = 0.79 cm for the conditions in Fig. |5J), the two skin 
layers near both walls are independent of each other, as 
the gap width is much larger than the nonlocality length 
I = Vt/u = 1.55 cm. As a result, for L > 5 cm the 
surface impedance does not depend on the gap width. 
In the opposite limit L <C 26, the electric field profile 
is linear E y = E y (Q)(l — 2x/L) and the plasma resistiv- 
ity is mostly determined by the first bounce resonance 
u> = fi;,. For uniform plasmas, only slow electrons con- 
tribute to the collisionlcss heating, because the resonant 
velocity corresponds to v x = uL/n = VrL/nl « Vt- If 
L <C 26 the number of resonant electrons and the total 
number of electrons in plasma are decreasing for smaller 
L, which leads to a smaller heating. For nonuniform plas- 
mas, the bounce resonance condition u> = Oj, can not be 
satisfied for any electron energy for small L < 2 cm and 
collisionless heating vanishes, see Fig- EI 

The maximum of plasma resistivity occurs at L ~ 26. 
Similarly to the case of short gaps, the real part of the 
surface impedance is mostly due to the first bounce reso- 
nance to — Qf,, and only slow electrons contribute to col- 
lisionless heating for uniform plasma, as v x ~ 28oj/tt = 
Vr26/irl < Vt- For nonuniform plasmas, the bounce 
frequency is higher than in uniform plasmas and most 
of electrons are in resonance. As a result the surface 
resonance plotted as a function of the gap width has a 
pronounced peak compared to the shallow maximum in 
uniform plasmas, see Fig. [21 

To compare the condition of the bounce resonance 
u) = Qf, with the corresponding condition of the tran- 
sit resonance u> = Vrk, additional simulations of the 
dependance of the plasma resistivity KeZ on the driv- 
ing frequency u> were performed for the fixed length 
L = 50 cm with given Maxwellian EEDF correspond- 
ing to the electron temperature T e = 5 eV, and vari- 
ous electron densities for a uniform plasma (without ac- 



counting for the electrostatic potential). In Fig.[3fa) the 
resistivity is shown as a function of the discharge fre- 
quency normalized by the inverse of the "transit" time 
of the electron pass through the skin layer Vr/<S m ax> 
where 6 max is the plasma penetration depth of the rf 
field for the maximal value of the resistivity defined as 
6 = f^ /2 dx\{E y (x)/E v (0)\. In Figure EJb) the plasma 
resistivity is plotted versus the driving frequency nor- 
malized by the thermal bounce frequency = Vrn/L. 
At high plasma densities (> 10 10 cm~ 3 ) the field pene- 
tration depth is much smaller than the discharge gap L 
and the effects of the finite discharge gap are unimpor- 
tant: the maximum of the surface impedance and, corre- 
spondingly, the most efficient collisionless electron heat- 
ing occurs at the "transit" resonance u) ~ Vr/6, as for 
the case of semi-infinite plasmas. At low electron densi- 
ties (and frequencies) (n e < 10 10 cm -3 ) the electric field 
penetration depth is of the order of the discharge dimen- 
sion 6 ~ L, and as a result the maximum of the surface 
impedance corresponds to the condition of the bounce 
resonance u> ~ ttVt/L. Note that the absolute maximum 
of the power dissipation occurs when both conditions for 
the bounce and transit resonances are met, which occurs 
for L = 7r<5 15]. This happens at n e < 10 9 cm -3 for the 
conditions of Fig. At high frequencies ui ^> irVr/6, 
the number of resonant electrons is exponentially small 
and collisionlcss heating vanishes. The resulting heating 
depends on the collision frequency and nonlocal effects, 
as described in Ref. []|. 

Accounting for the ambipolar potential enhances con- 
siderably the plasma resistivity for the aforementioned 
case. As shown in Fig. the pronounced maxima of 
Reif appear for the frequencies that correspond to inte- 
ger multiples of the bounce frequency u> = nil^x, because 
accounting for the electrostatic potential yields a larger 
number of the resonant electrons. 

Self-consistent calculations 

To investigate the behavior of discharge parameters 
under the condition of a bounce resonance, the full self- 
consistent simulations of the EEDF, rf electric field, and 
ambipolar potential for fixed surface current have been 
performed for 13.56 MHz rf driving frequency. Figure 
shows the dependence of the resistivity of the dis- 
charge plasma on the discharge dimension. The calcula- 
tions have been performed for discharge gaps from 3 cm 
to 10 cm (the discharge can not be sustained for gaps 
smaller then 3 cm at a given pressure of 3 mTorr). It 
can be clearly seen that the resistivity of the plasma 
sharply increases for the discharge gap corresponding to 
the bounce resonance (about 3 cm). The self-consistent 
electrostatic potential and ion-electron density profiles 
are plotted in Fig.|§fa) for two different discharge lengths 
- 3 cm, corresponding to the bounce resonance condition, 
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FIG. 4: The plasma resistivity as a function of the driving 
frequency normalized by the thermal bounce frequency for 
different electron densities with and without the ambipolar 
potential <p(x) — 5x (2x/L— l) 4 eV and for the same discharge 
parameters as in Fig. 3. 



and 10 cm, corresponding to the non-resonant width. 
These graphs show that the electron density at the center 
of the discharge is larger for the 3 cm resonant gap than 
for the 10 cm non-resonant gap. Note that if the power 
transfer efficiency, or the surface impedance, were the 
same, then the total power transferred into the plasma 
would also be the same and the plasma densities would 
be equal, due to energy balance. In our case the surface 
impedance for 3 cm gap is considerably higher, what cor- 
responds to the higher plasma density. 

The electron energy distribution function and the dif- 
fusion coefficient in energy space are shown in Fig. EJb). 
Figure HHb) shows that the energy diffusion coefficient 
is larger for the 3 cm gap than for the 10 cm gap for 
electron energy less than 15 eV. This results in more 
effective electron heating, leading to the larger plasma 
resistivity shown in Fig. [5] The steady-state electron 
energy distribution function is governed by the follow- 
ing processes: the collisionless electron heating in the rf 
electric field, inelastic collisions with neutrals, and re- 
distribution of energy among plasma electrons due to 
electron-electron collisions. We see in Fig. E{b) that the 
EEDF shape is similar to the two-temperature EEDF 
4] with the temperature of the tail of the distribution 
being lower than the temperature of the main body of 
the EEDF, corresponding to the onset of inelastic col- 
lisional losses. For a 3 cm gap, corresponding to the 
bounce resonance condition, the electron temperature of 
low-energy electrons (less than the excitation potential 
11.5 eV) is much higher than for the 10 cm non-resonant 
gap. This effect is similar to the plateau formation on 
the EEDF governed by collisionless heating in the finite 




FIG. 5: The results of self-consistent simulations for the 
plasma resistivity and the electron temperature (defined as 
2/3 of the average electron energy) at the center of the dis- 
charge as functions of the size of the discharge gap and 
for a given surface current / = 5 A/m, rf field frequency 
Lu/2-K — 13.56 MHz and argon pressure P — 3 mTorr. 



range of electron energies Under the conditions of 
Fig. El this plateau is not well pronounced, because it is 
smeared out by electron-electron collisions. Additional 
simulations have been performed also for the discharge 
frequency 6.78 MHz, which is half the driving rf field 
frequency considered above. For the lower driving fre- 
quency the first bounce resonance shifts in the region of 
larger L, as shown in Fig. [7| Figure [7\ shows the calcu- 
lated resistivity for two different surface currents, 1 A/cm 
and 5 A/cm. One can see that the positions of the re- 
sistivity maxima corresponding to different surface cur- 
rents are shifted relatively to each other. The larger sur- 
face current corresponds to a larger power transfer into 
the plasma according to Ea. H37|l and results in a higher 



plasma density (n 



= 2x 10 10 and n e = 7x 10 11 , 



respectively). The higher discharge plasma density, in 
turn, leads to a smaller skin depth. Correspondingly, the 
position of resistivity maximum shifts into the region of 
smaller discharge gaps of the order of the skin depth. The 
electron energy distribution functions for 6.78 MHz are 
plotted in Fig. [5] for the surface current 1 A/m and for 
two different lengths, resonant 9 cm and non-resonant 
15 cm. The phenomenon of plateau- formation on the 
EEDF is clearly seen for the bounce resonance condition 
for L = 9 cm. 



CONCLUSION 

The analysis of the properties of inductively cou- 
pled discharges clearly shows the phenomenon of the 
bounce resonance. Self-consistent simulations of the dis- 
charge resistivity and electron energy distribution func- 
tions demonstrate the significance of the explicit account- 
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FIG. 6: The results of self-consistent simulations for the dis- 
charge gap L — 3 cm corresponding to the bounce resonance, 
and L = 10 cm for the same conditions as in Fig.5 (a) the elec- 
tron density and ambipolar potential profiles, (b)the electron 
energy distribution function EEDF and the energy diffusion 
coefficient D E (e) profiles. 



ing for the non-uniform plasma density profile and the 
correct form of the ambipolar electrostatic potential. En- 
hanced electron heating and larger plasma densities (for 
a given current in the coil) can be achieved if the low- 
pressure ICP discharges are operated under the condi- 
tions of the bounce resonance. 
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